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I. INTRODUCTION 



A long-standing problem in the theory of polymer molecules in solution is the effect 
of intramolecular forces on the shape and size of an isolated polymer. The intramolecular 
forces are usually assumed to be of van der Waals' type, consisting of strong, short- 
range repulsive and weak, long-range attractive interactions. Depending on temperature or 
solvent composition, the macromolecule forms either extended structures or collapses into a 
dense globule 0. The transition between these two states, which is called 6 transition, has 
attracted a great deal of attention, both theoretically [|8Hl3| and experimentally partly 
because of its connection with protein folding |[T^-[16|. 

For the bidimensional linear polymers many theoretical tools, such as conformal invari- 



ance [17,18| and Coulomb gas [19| methods, can be successfully employed and considerable 
progress has been made in understanding the nature and the critical exponents at the col- 



lapse transition |20,21 



In contrast much less is known about branched polymers. In lattice statistical mechan- 
ics they are described by lattice animals, graphs of connected occupied sites on a lattice. 
Some examples are shown in Fig. 1. Contact interactions are introduced between nearest 
neighbour sites which are not immediately linked by a bond and solvent interactions are in- 
troduced between an occupied site and any nearest neighbour empty sites. In d dimensions 
lattice animal exponents are related to the exponents of the Lee- Yang edge singularity of the 
Ising model 'm.d—2 dimensions |]22|. However, this elegant theoretical approach can give only 
the bulk entropic exponent for the extended phase. Moreover, branched polymers in two 



dimensions are in general not conformally invariant [Q. This model can be reformulated 



as the g — > 1 limit of an appropriate g-state Potts model [0,^]. However this formulation 
has not given any exact results but only allowed the conjecture that lattice animals 
could present two distinct branches of 9 behaviour, separated by a percolation critical point 



While there is little exact or analytical information available for this problem, a large 
amount of numerical work has been done using Monte-Carlo |]27|-p9|, exact enumeration 
||30|,|8|j9|l and transfer matrix methods. The results of these analyses are controversial, 

however, and it is this which triggered the work to be described here. 
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In order to describe the purposes and the goals of this paper let us first introduce some 
notation. 

Let n, b, s, k denote the number of sites, bonds, solvent interactions and contact interac- 
tions of the lattice animal. We are working on a square lattice with coordination number 
7 = 4. Then the relation^ 'jn = 2b + 2k + s implies that three of the four variables are needed 
to specify the problem completely. The statistical mechanics of the model is most conve- 
niently described by a generating function. This has recently been defined in the literature 
PJTTH in terms of two sets of variables and it is useful to give both of them here. 

n n,b,k 

= e'^°"Z„(/5i, (32) = E <'n,s,k exp {(3on + (3,s + (32k) (1.1) 

n,s,k 

where an,b,k stands for the number of graphs with given n, 6, k and x, ?/, r and /3o, A, (32 are 
fugacities. The relationship between the two sets of variables is given by 

a; = exp(/5o + 7A) , y = exp (-2/5i) , r = exp (/Js - 2/5i) (1.2) 

The infinite lattice animal can be described either by taking n — oo or by fixing one of the 
three fugacities in terms of the other two 0, say x = Xc{y,T) or Po = /3oc(/?i, /?2) (see also 
below) . 

The compact and extended phases can be distinguished from the scaling behaviour of 
the mean radius of gyration, (Rn), with the number of monomers A^, viz. 

{Rn) ~ iV^ (1.3) 

The phase diagram is sketched in Fig. 2. First, there is a (second-order) transition from an 



extended state (with u ~ 0.6408 in two dimensions ^2[) where the different branches of the 
lattice animal are widely separated, to a compact state {u = 1/d). The locations of this 
extended-compact transition which follow either from exact enumeration P,|33| or transfer 



matrix calculations p2|JTT[| are in good agreement with each other. However, there is an 
important discrepancy, which is the subject of this paper. While the enumeration studies 



^In two dimensions, this is valid for 7 < 4. 
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33[] suggest the existence of an additional transition between two distinct compact phases, 



no sign of such a transition was found using the transfer matrix |]rT] . 

The reason why a compact-compact transition could in principle be present can be better 
understood by looking at Fig. 1. Among the compact configurations (a,b,c) which can fill all 
space, there are two different types. The first is a contact-rich (and almost linear) branched 
polymer (Fig. la,b) and the second type is a cycle- rich branched polymer (Fig. Ic). Contact- 
rich configurations are expected to provide the dominant contributions to Z in the r — oo 
region, whereas cycle-rich configurations are expected to dominate as ?/ ^ oo. The question 
is whether or not there exists a sharp transition between the compact configurations.^] 

Our aim in this paper is to further investigate on the existence or otherwise of this 
transition, by considering some approximations which render the model analytically solvable. 
This will be done using an exact mapping onto an extended g-state Potts model, which will 
be defined in the next section. The results of the approximate analysis of the extended 
g-state Potts model will then be reinterpreted in terms of the lattice animal. Some of these 
results were already stated without derivation in [Q. In section 3, we present the exact 
solution of the extended Potts model on the Bethe lattice and derive the phase diagram. 
This calculation will be supplemented in section 4 by an analysis of the fixed-point structure 
of the model resulting from the Migdal-Kadanoff renormalization group equations. In both 
schemes our results support a phase diagram with a single extended and a single compact 
phase. These are separated by a line of 6 transitions, divided into two sections which fall into 
different universality classes and which are separated by a higher order multicritical point, 
which coincides with critical percolation. In section 5, we illustrate further the role of the 
competition of the contact-rich and cycle-rich states by studying two-dimensional directed 
branched polymers, where the cycles are absent. We find a single line of 6 transitions in the 
directed percolation universality class. Finally, in section 6 we present our conclusions. 



similar kind of transition has recently been discussed in the polymer literature |12,34| and 
it displays some analogies with phenomena occurring in polymer gels, for which a multiplicity of 
phases can be realized in the collapse process 



II. EXTENDED POTTS MODEL FORMULATION 



We are interested in the phase diagram of infinite branched polymers which are defined 
by the generating function Z in eqn (|1 . 1|) . Two existing results can serve as a guide. Firstly, 
in the case of strong embedding (no contacts allowed), the extended-compact transition is at 
y ~ 6.48 and r = [^]. Secondly, it was shown [0,[TT]| that along a certain line in the x, y, r 
(or /5o, /?!, [32) space the branched polymer can be mapped onto percolation. The percolation 
critical point corresponds to the following point on the 9 transition line 

yp = 2 , Tp = 2 (2.1) 

(or /?i = -iln2,/52 = 0). 

Our approach is based on the known |2^j37|j2^j25[] exact relation of the lattice animal 



with the extended g— state Potts model with the classical Hamiltonian 

n = -jY. - ^ E ^^.^^^.i -hY. ^^.a (2-2) 

where cxj = 1, 2, . . . , g and the fugacities are expressed in terms of J,L,H as 

X = exp[-H - -f{J + L)] , y= (e^ - l) e*^+^ , r = 6-^+^ (2.3) 

The relation of this extended Potts model with the lattice animal problem is 

Z = /potts = hm I- In Z , Z = Y,e''^ (2.4) 
9-1 dq ^ 



where /potts is the Potts model free energy. The connection in eq. ( |2.4| ) can be established 
by considering the graphs in the high-temperature expansion of the Hamiltonian ( p.2|) and 
matching the weights in such a way as to reproduce the branched polymer generating func- 
tion. Thus, calculating the phase diagram of the Potts model gives information about the 
phases of the animal problem. This correspondence is exact. However, approximations are 
needed to treat the extended Potts model. Having solved the extended Potts model in the 
approximations stated below, we shall reinterpret the results for the extended Potts model 
phase diagram in terms of the lattice animal collapse transition(s). We shall pursue two 
approaches. 
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1. The exact solution on the Bethe lattice should give the correct topology of the phase 
diagram for sufficiently large spatial dimensions d. While the Bethe lattice solution is 
expected to become exact in the d ^ oo limit and should have the same universality 
properties as a mean field treatment (see |^8[ for a detailed discussion), it has the 



advantage over mean field theories that the self-consistency equations to be dealt with 
are considerably simpler. 

We point out that for the original branched polymer problem, a Bethe lattice solution 
is of little interest, since contact interactions are eliminated on the Bethe lattice by 
construction. Contact interactions may be kept by considering the lattice animal on 
a cactus lattice, which is the dual lattice of the Bethe lattice, see On the other 



hand, contact interactions are not projected out by our treatment of the Bethe lattice 
approximation in the extended Potts model, as we shall see below. 

2. We supplement this by an analysis of the fixed point structure of the model, using the 
Migdal-Kadanoff approximation, which should become reliable in d = 1 + e dimensions. 



as reviewed in Our finding that the structure of the phase diagram is the same 
in both approximations makes it plausible that the topology is valid for all dimensions. 
This result is the more remarkable since it is known, for example in three-component 



lattice gases [^, that the phase diagrams found from mean field theory and Migdal- 



Kadanoff approximations may in general be distinct. 

Our aim is to find the phase diagram of infinite lattice animals through a study of the 
equivalent extended Potts model. It is well known that taking the limit n — oo corresponds 



to finding critical points of the Hamiltonian (|2.2|) . Transitions between different infinite 
polymer states then correspond to multicritical points within the critical manifold of the 
extended Potts model, see 101. 



We now recall a relationship |TT| which will later be used to simplify the calculations on 
the Bethe lattice. In principle, we are interested in the lattice animal free energy F. In the 
canonical ensemble (with n fixed), it is rigorously known that 



F(/3i,/?2) = Jim n-MnZ„(/?i,/52) (2.5) 
exists and is a convex function in (3i and j32- We now apply eq. ( p.5|) to fix x in the grand 



canonical ensemble such that we describe the infinite lattice animal. Using eqn ( p.5|) , the 
animal generating function takes the asymptotic form Z ~ Z^nl^^^^)"; which diverges at the 
critical point Xc given by 

F=-Inxc. (2.6) 

Thus, working in the grand canonical ensemble, it is sufficient to find as a function xj^y^ r) 
of the other fugacities y and r. Once the critical value of x is obtained as a function Xc(y, t), 
the canonical free energy follows from eqn ( ^^.61 ). 

III. BETHE LATTICE SOLUTION 

We now describe in detail the exact solution of the g — 1 extended Potts model on the 
Bethe lattice. 



A. Generalities 



Consider a Cayley tree (Figure 3) with 7 = 3 nearest neighbors for each site. Following 



the terminology of |38|, by Bethe lattice solution we mean the behaviour deep inside the 
Cayley tree. Let cxo denote the central spin of the Bethe lattice. Then the extended Potts 
model partition function is 

^ = ^e^ViEng„(a„|,0)) (3.1) 

o-o {s\ j=l 

where s^^^ denotes the spins in the j*^ sub-branch and 

Qn (O-Q \s) = exp {JSaQ,si + L6ao,lSs^,l + Hd^g^i) 

■ exp ( JEX,.. + ^EX,i5s,,i + i^EX.i ) (3.2) 
where J2' is a sum over one sub-branch of which the first site is si. 



B. Derivation of the self-consistency equations 



We define, following Baxter |^ 
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9n{(^o) =Yl'Qr^i^o\s) (3.3) 

{s} 

where n refers to the number of iterations performed in the construction of the Cayley tree. 
Recursion relations for the Qn are 

<? 

9n{cro) = J2 ^^P (■J^'^o,si + LS^o^iSs^^i + H5s^^i) ign-iisi)f . (3.4) 



Having solved eqs. (|3.4| ) and thus obtained the gn, the thermodynamics is determined com- 
pletely. 

The analysis of the eqs. (|3.4| ) can be simplified, in analogy to the procedure described in 
4^ , by introducing the variables 



„ _ gnjfTo ^ 1,2,3) _ gn{2) _ ^n(3) , . 

^n(l) ^n(l) 9nO-) 

which, as will become apparent later, are sufficient to describe the full thermodynamics. 
Our choice of variables was motivated by the mean field treatment of the ordinary g-states 



Potts model in an external symmetry-breaking field [42]. There it was shown that for the 
distinct ground states corresponding to different phases the order parameter either has the 
same value for all its q components or that at most the order parameter component which 
is directly coupled to the field may be different from the others. In ( p.5| ) we generalize 
that result in order to allow for the possibility of a compact-compact transition besides the 
expected 9 transition. As n oo, the variables in ( p.5|) tend to fixed point values S, T, Z, 
which can be determined from the recursion relations ( p.4[ ). At this stage, we take the g — >■ 1 
limit (see (|2.4|) ) and find the self-consistency relations 



-1-^-3 



X T 



+ {y/T + l)T' + Z^-2E^ 



a;-ir-2 + T2 Z2 - 2S2 ' ^ ' ^ 

x~W-^ + + jy/r + 1)Z^ - 2'E? 
~ a;-ir-2 + T2 + ^2 - 2~? ' 

The equations ( |3.6D can be decoupled after some transformations of the variables. We then 
find that there are four cases which must be considered separately. As we now show, in each 
case the problem of solving a set of three coupled equations can be reduced to the solution 
of a single quartic equation. 
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1. Case A: S = T = Z. 

The fixed point equation is a quadratic equation in S, with the solutions 



„ 1 ± yi - Axy 

-± - 7, I'^.rj 

zxyr 

From eqn (|3.6|), it is clear that S_ is the stable solution which, using eqn ( p.4| ), deter- 
mines the lattice animal generating function 

Z = xT^El (3.8) 

These solutions are real provided Axy < 1. Criticality corresponding to an infinite 
lattice animal is recovered at 

X. = 1 (3.9) 

2. Case B: H = Z ^ T. 

We introduce the variables 

u = E + T , v = T-E. (3.10) 

For V = 0, case A is recovered. For v ^ 0, & factor v cancels from the fixed point 
equations (|3.6|) giving 



v = --, (3.11) 

T XT"^ U 

x^yr^ {u^ - 4m^) + xt^ (4(r - 1) - xy^) + Ixy^Tu - y = 0. (3.12) 



3. Case C: S ^ T = Z. 

We introduce again u and w defined by (|3.1CI|) . For t> = 0, case A is again recovered. 
For f 7^ 0, a factor v cancels giving 

v = — -^-, (3.13) 
2r 2xT^u ^ ' 



x"^ 



yr^ (u^ - 4u^) + sr^ U{t - 1) - -xy^\ + ia;yVu - - = 0. (3.14) 
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4. Case D: S ^ T ^ Z. 

We define new variables 



a = T + Z , b = T-Z (3.15) 



and rewrite the fixed point equations ( |3.6|) in terms of E,a,b. 6 = corresponds to 
case C. For 6 7^ 0, a factor b cancels giving 

62 = -a2 + 4S2 + ^a--^ (3.16) 



together with 



T XT^ 



^a' = —^- A-E? + {y/r + 2) f + 22^ - 

r XT-^ \T XT'' 

y ~- L + ^■^a ^ y_ 

-3 ^ ^ rrn-'^ 



-ar. = — r- + -^^ + -a -. (3.17) 



Provided S 7^ 1, we solve the second of eqs. ( |3.17|) for 



^(S-l)a = ^H2 + ^ (3.18) 

and find 

x^y V (S^ + {y/T - 2)E') - xy^T^l + xi/r^)^^ 

+ x?/r(?/(l + r) + 2r(l - r))H - (1 - r)^ - xyV = (3.19) 
Note that the relation between S and Z remains unspecified. 



C. Analysis of the self-consistency equations 



The problem has been reduced to the solution of three quartic equations ( p.l2|) , (p.l4|) 



and (|3.19| ). Further analysis is greatly simplified by numerically checking that all solutions 



in Case D where S, T and Z are a priori different, reduce to one of case B, with S = Z. 
As we shall show below, this implies that there are only two distinct phases for the infinite 
lattice animal on the Bethe lattice. In particular, that means if any two branches of the 
lattice animal generating function calculated in the cases A,B,C meet, all three solutions 
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have to coincide. It follows that the topology of the phase diagram calculated in the whole 
(S, T, Z)-space is the same as would have been found when only considering (S, T)-space, 
with T = Z. We expect the same reduction to occur in even more generic setups than done 
in eq. (p.5|). 

We only have to consider the cases A,B,C. We introduce a new variable 

p = UT , u = E + T. (3.20) 



We recall from (|2.6|) that in order to get the infinite lattice animal free energy in the canonical 
ensemble, it is enough to know x as a function of y and r. Solving for x = Xc{p', y, t) rather 
than solving for p (or u) has the further advantage of going from quartic to quadratic 
equations, p now plays the role of an order parameter and will be fixed by maximizing x{p) 
with respect to p. We find 

1. Case A. We have 

Aip- y,T):= x(^) ip; y, r) = ^^^^ (3.21) 

yp 

2. Case B. The equation is 

x^iyp^ - Ayrp^) + x (4(r - l)p^ - V + 2y'^p) -y = (3.22) 

with the solutions 

Bi{p;y,T):=x^^\p;y,T) (3.23) 

,2^ _ Ao,2 



2(1 -T)p-y^ ^ ^p(4r2-8r + 4 + y2)_4y2 



3. Case C. The equation is 



x\yp^ - Ayrp') + x ( 4(r - l)p' - -xyY + -y'qj - | = (3.24) 



with the solutions 



C±{p;y,T) ■.= x^f\p;y,T) = ^B±ip;y/2,T). (3.25) 
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To get the equilibrium physics, we have to minimize the lattice animal free energy F{y, r) 
(see eqs. ( P3| , ^.6D ) with respect to p or equivalently to maximize x{p;y,T). For case A, we 
have 

'2(p-2)\ dF p-A d'^F 



In X = — In 



p'^y I ' dp p{p — 2) ' dp'^ 



= l>0 (3.26) 

p=4 ^ 



and F has a single minimum at p = 4. Indeed, F = ln{4y) evaluated at p = 4 is concave in 

y- 

It remains to be seen under which conditions this local minimum of F^^^ is the absolute 
minimum^ of F. /^From the above it is clear that if any two of the solutions A,B,C meet, 
we have in fact a simultaneous meeting of all three of them. This common meeting point 
occurs at p = 4 

A{4;y,T)=B44;y,T)=C4A;y,T) = ^ ; r > 1, 

Ay 

A{4; y, r) = B+{4; y, r) = C+(4; r) = ^ ; r < 1. (3.27) 
In addition, we have checked that 



C.{p-,y,r)>B_{p;y,T). (3.28) 

For an example illustrating the solutions A, B^,C^, see Figure 5. This implies that there 
are just two distinct phases of the lattice animal which are described by the two solutions 
A{p;y,T) and C-{p;y,T). Furthermore, we want solutions such that the fugacity x is real 
and positive. This yields the condition 

Case A 

(3.29) 

i^^^^ Case C 




(r-l)2 + (y/4)2 

On the other hand, the unstable solution S+ from case A always corresponds to p > 4. 
Furthermore, we argue in the appendix that for p > 4 the solution C_ is not a stable 
solution of the self-consistency relations (|3.6| ). We therefore must have 



^It is enough to consider B- and C_ , since the other two solutions , (7+ have no maxima with 
respect to p. 

12 



p < 4 (3.30) 

To get the transition lines between the two phases, note that dC-/dp = if and only if 
r = 2. Furthermore 



C_(p; y, 2) has at p = 4 a 



maximum if y > 8 

turning point if ?/ = 8 (3.31) 
minimum if y < 8 

To locate the transitions, one must distinguish two cases. First, for y < 8, one has a first 
order transition. It is given by the conditions 

C^{p;y,T) = A{4;y,T) = 

Ay 

dC 

--{p;y,T) = (3.32) 



dp 



< 



Qp2 

At this transition point, p jumps from its value p = 4 for r small to a new value Pc{y, t) < 4. 
The location of the line and Pc have to be found numerically. 

Second, for y > 8 and r > 2, numerical studies show that C^{p; y, r) has a maximum at 
some value p' < 4. For all p < 4, C_ > A and the state described by p' is the equilibrium 
ground state. As p ^ 4, the two solutions approach each other continuously. 

Finally, for y = 8 and t = 2, the second order line ends in a tricritical point. The two 
transitions are shown in Fig. 6. We also show the jump in p, Ap = 4 — pdy, t) and observe 
that Ap ~ (8 — y) for y 8^. 



D. Summary 

Studying the phase diagram of the extended Potts model in the g — >■ 1 limit on the 
Bethe lattice, we find that there are two distinct phases (A and C above) on the manifold 
corresponding to the infinite lattice animal. The transition between them is first order for 
y < 8 and second order for y > 8. There is a tricritical point at ?/ = 8, r = 2. The Bethe 



lattice calculation corresponds to a mean-field calculation |38|. Detailed comparison |38 



of the exact solution of the simple q— states Potts model (that is, L = H = 0) in two 
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dimensions with the Bethe lattice solution suggests that the latter faithfully represents the 
model behaviour for infinite dimensionality. We expect the same for the more general models 
treated here. 



IV. MIGDAL-KADANOFF RENORMALIZATION GROUP 



In the last section, we studied the extended Potts model ( |2.2| ) using an approximation 
which should become more reliable with increasing dimensionalty. We now wish to com- 
plement this calculation by a Migdal-Kadanoff renormalization group study. The Migdal- 
Kadanoff approximation preserves the self-duality property of several two-dimensional mod- 
els (in particular the g— state Potts model, see [^0|,^) and thus yields the phase diagram 



exactly. The bond-moving approximation involved is exact at zero temperature and its pre- 
dictions for exponents are exact to order e in d = 1 + e dimensions, as reviewed in detail 
in 



40[| . We therefore hope that the renormalization group calculation will provide reliable 
information on the fixed point structure for low dimensions. 

The Hamiltonian ( p.2|) has already been investigated using a Migdal-Kadanoff renormal- 
ization group by Coniglio He used a rescaling factor 6 = 2 in dimensionality d = 2 
and considered q = 1, rather than q 1^. He identified four non-trivial fixed points which 
correspond to (1) an extended lattice animal phase, (2) a compact phase, (3) a percolation 
point and (4) a tricritical point analogous to a point. These are not sufficient to describe 
the topology of Fig. 2 and therefore it seemed useful to look again at the calculation with 
particular emphasis on the possibility of the occurrence of a compact-compact transition. 

Here, we shall consider the case d = 1 + e. Furthermore, as will become apparent, the 
limit g — s> 1 has to be taken with care, because the fixed point structure changes at g = 1. 
For d = 1 + e and 6 = 2 the Migdal-Kadanoff recursion equations for the extended Potts 
Hamiltonian ( p.2|) are 

. fi(l±ft(iziMV' (4.1) 
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where 



^ = exp(i7/2 - J) , T] = exp{-J) , p = exp{L + H) (4.4) 



Eqs. (4J)-(13) were obtained by performing a one-dimensional decimation followed by 



bond- moving PO| , p | 



The fixed point structure that follows from the recursion equations ( [4.1[ - [4.3|) is compli- 
cated and g-dependent. However, using the q = 2,3 cases as a guide in the two limits 
of interest to us (g ^ 1 and e ^ 0) a clear pattern appears and 14 fixed points can be iden- 
tified. These are listed in Table ||, together with the eigenvalues of the Jacobian evaluated 
at the fixed points. We also illustrate the fixed point structure in Fig. 6. 

Note that that the fixed points naturally divide into three groups. The first group 
contains the fixed points B, C and D for which p* = cxd. The second group, which is the 
one we shall be interested in, contains the fixed points E,F,G,H and 1 which all have p* ~ 1. 
Finally, the third group contains the fixed points with p* ^ 0, namely A and the cluster of 
fixed points around y — 1, J,K,L,M and N. For q 1, the last five points merge to a single 
fixed point. 

To understand the physics, we begin by singling out the completely attractive (or trivial) 
fixed points. Reading from table |I| the fixed points B,C,D and A,J,N only have eigenvalues 
which are smaller than one. Hence they cannot describe the critical behaviour of the infinite 
branched polymer. Next, we consider the cluster of five fixed points near to (^*,?7*,p*) ~ 
(0,1,0). One (M) has two relevant eigenvalues, whilst two (K and L) have one relevant 
eigenvalue, with an associated exponent y = d. From the eigenvectors, we have checked that 
the renormalization group flow from the vicinity of one of the fixed points M,K,L always 
goes towards the trivial stable fixed points J and N. Since these five fixed are going to merge 
into a single fixed point in the q —>■ 1 limit and none of them has three relevant eigenvalues 
as would be required from a percolation fixed point, we conclude that they cannot govern 
any part of the phase diagram of the branched polymer problem. We thus find that the fixed 
point structure depends on carefully taking the g ^ 1^ limit, rather than simply setting 
g = l. 

It remains to consider those fixed points which are close to the p* = 1 plane. Two of 
these (E and I) are independent of both q and e and have one relevant eigenvalue. For the 
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fixed point E, the relevant direction is characterized by the exponent l/u = d = 1 + e. This 
is characteristic of a compact object and we therefore conclude that this fixed point controls 
the compact phase of the infinite branched polymer. For the fixed point I, the single relevant 
eigenvalue is l/u = d — 1. Since u > 1/d, we expect this to describe a non-compact phase 
and the fixed point I to govern the extended phase of the polymer. 

The fixed points F, G and H are dependent on q and e, but for g < 2 all appear in the 
plane p = 1, to leading order in e. As d decreases they move towards the point E. The fixed 
point H has three relevant eigenvalues and therefore represents the percolation fixed point 
(which is realized for H = and L = in ( |2.2| ) and ioi d = 2 the Migdal-Kadanoff approach 
is known |40,4^ to correctly reproduce the exact critical point). F and G are tricritical 
points and govern the renormalization group flow along the two critical lines leaving the 
percolation point H. We point out that the fixed point G is only found when the limit 
g — s> 1+ is carefully taken. However, if one simply puts q = 1, that fixed point is missed, in 
agreement with the earlier work of Coniglio ||26|| . 

Thus we predict a line of 6 points which contains two segments falling into two distinct 
universality classes and meeting at the percolation point. The 6 line separates the extended 
phase from a single compact phase. 



V. DIRECTED LATTICE ANIMALS 



So far, in this paper we have considered the structure of the phase diagram of the infinite 
branched polymer. Although the interplay of the different ground states (see Figure 1) in 
the contact-rich and contact-poor region of the compact phase might at first sight suggest 
the existence of a compact-compact transition, this was not supported by our (albeit ap- 
proximate) calculations. To provide some insight into the role of the different ground states, 
we now wish to consider a model where the spiralling or folding ground states shown in 
Fig. la,b are eliminated. 

For this reason, we now consider directed branched polymers. This model is defined in 
complete analogy with the system discussed so far, but requiring that the bonds, solvents 
and contacts are directed. We use a square lattice with the preferred direction along a lattice 
diagonal. Animal configurations are allowed only if (1) they start from a single site (2) all 
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bonds have a component in a preferred direction and (3) the solvents and contacts counted 
in the generating function have a component in the preferred direction. The special case 
where contacts are forbidden (that is, the limit /?2 ~^ — oo) was studied previously by Dhar 
44 . 

The interest in this model arises because due to the directedness condition ground states 
of the undirected lattice animals in the contact-rich region — — oo, — > oo) are not 
allowed. Indeed, the phase diagram of the directed model (Fig. 7) is quite different, as we 
shall now show. 

For directed lattice animals, it is always possible to reduce the problem to the "strong- 
embedding" {P2 —00) case. To see this, consider a strongly-embedded lattice animal 
(where all neighboring sites are connected by bonds). Strong embedding is broken when 
some of the bonds are replaced by contacts. We now ask how often these replacements of 
bonds by contacts can be made. 

Because of the directedness, it is clear that bonds can be replaced by contacts only if the 
strongly-embedded lattice animal contains closed loops. Furthermore, each loop contains 
exactly one "final" site. Only one of the two bonds leading to this final site may be replaced 
by a contact because of condition (1). Now, consider a directed lattice animal with c loops 
and denote by g{c, k) the number of ways of replacing exactly k of the bonds by contacts. 
Obviously, g{c, k) = ii k > c, since no loop can be cut twice. It is also clear by inspection 
that g{c, 0) = 1 and g{c, 1) = 2c. Also, for > 2 the recursion formula 

g{c, k) = 2g{c -l,k-l)+g{c- 1, k) (5.1) 

is valid. To see this, imagine placing the first contact onto the directed lattice animal. Single 
out one of the loops. Either this contact cuts this loop, which can be done in exactly two 
ways, and one then has to place k — 1 contacts onto an animal with c — 1 loops or this loop 
is not cut at all and all k contacts are placed onto the remaining c — 1 loops. Defining the 
function 

fic;P) = i:gic,k)e^' (5.2) 

A:=0 

it is easy using (|5.1j ) to show that for all c > 

/(c; (3) = {l + 2e^) /(c - 1; /?) = (l + 2e^)' (5.3) 
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Consider now the directed lattice animal generating function. On the square lattice, 2n 
b + k + s. Using the Euler relation c = b — n + 1 and 



g 

= J2a{n, 0, s)e*"e'5i^ (l + 26^^)' (5.4) 
g 

where a{n, k, s) is the number of directed lattice animals with n sites, k contacts and s 
solvents. Q runs over all directed lattice animals and Q over all strongly embedded directed 
lattice animals. In the strong-embedding case, k = 0. Thus c = n — s + 1 and 



Z={l + 2e^^) Y.a{n, 0, s) (e* (l + 2e^^))" (- 

n V 



+ 26/32 



(5.5) 



g 

This is the strongly embedded directed lattice animal, with the effective solvent fugacity 

= -——r (5.6) 

1 + 2e^2 ^ ^ 



The strong-embedding model can be mapped onto directed percolation WM. We have thus 



shown that the fully directed lattice animal has a single extended-compact transition which 
is in the directed percolation universality clas^. In contrast to the undirected model, there 
is no change of the values of the critical exponents along the transition line. 



VI. CONCLUSIONS 

In this paper, we have investigated the phase diagram of the extended Potts model ( p.2|) 
in the g ^ 1 limit. We have concentrated on the structure of the manifold of critical 
points, because this is related to the phase diagram of an infinite branched polymer. Using 
Bethe lattice and Migdal-Kadanoff renormalization group calculations, we find a line of 
6 transitions between the extended and the compact phase. The 6 line consists of two 
segments which are controlled by different fixed points. Their meeting point coincides with 
the percolation critical point. Although these results were obtained from rather drastic 



^We have also checked that the relationship with directed percolation can be established directly 
from the fugacities given to the elementary processes of the directed animal model. 
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approximations, we expect the results coming from these to be rehable for large (Bethe 
lattice, d —* oo) or small (Migdal-Kadanoff, d = 1+e) dimensions, respectively. Furthermore, 
since the qualitative structure of the phase diagram is the same in both cases, it is plausible 
that the topology of the phase diagram is indeed independent of the dimensionality. Our 
results on the structure of the 6 line agree with previous numerical calculations using either 
graph enumeration or transfer matrix techniques. 

Our calculations give no evidence for a further transition between two compact phases. 
While that conclusion does agree with numerical transfer matrix calculations, it is in dis- 
agreement with analysis based on graph enumeration. We thus expect the crossing-over 
between the contact-rich ground states Figure lab and the bond-rich ground state Figure Ic 
to be smooth. While intuitively a phase transition between these ground states might ap- 
pear plausible at least in two dimensions because a ground state as the one in Figure la 
could be seen as a vortex while the state in Figure Ic is not, the nature of that transition 
is less apparent in three and higher dimensions. However, graph enumeration predicts ||3^ 
a compact-compact transition in three and more dimensions as well. 

This does not mean, however, that the different ground states of the model in different 
regions of the compact phase have no physical role. This is clearly illustrated by the case of 
directed branched polymers, where the spiralling states (Figure lab) are absent. Then the 
multicritical point on the 6 line found for the undirected branched polymer disappears and 
the entire 6 line is in the directed percolation universality class. 
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APPENDIX 



Here we study the stability of the solutions of the self-consistency equations ( |3.6| ). These 
arise from recursions of the form 

Xn+l = f{Xn) (A.l) 

and it is well known that a fixed point solution x* = lim^^oo Xn is stable under small 
perturbations if and only if |/'(a;*)| < 1. Consider first case A. Then we have the recursion 

S„+i = /(S„) = - + xi/rS^ (A.2) 
r 

With the solutions E± from eq. (|3.7|) we have 



/'(H±) = l±^l-4xy (A.3) 

and thus S_ is the stable and S+ the unstable fixed point. Since p± = 2rS±, it follows that 
if r ~ 2 as is the case for the A/C transition, we have p > 4 for the unstable case. Next, we 
look at case C. We introduce u, v from ( p.lOj ) and are interested in the case when v is small, 
that is, near to the 6 transition. Then the following fixed point equation for u holds 

2 1 

■u ~ — h -xyru^ + 0{uv) (A. 4) 

r 2 



which the same relation as for 2S. From (|3.7| ) it then follows that if m > m+ > l/(x?/r) 



it is an unstable solution. On the other hand, in the vicinity of the A/C transition we 
have from (|3.9| ) that at least in the transition region, all solutions which correspond to 
p = Tu > 1/ {xy) ~ 4, will be unstable. They are eliminated by requiring that p < 4. 
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TABLE I. Fixed points of the Migdal-Kadanoff recursion relations. The following ab- 



breviations are used: a = ^exp(— 1/e), (3 = exp ^— ^2^2 ) • 



\q-l)\u{0)-2H2-q) 



^ ln(2-9) . _ 
2-q ^ q-2 

2% Ai = (izHlMHzi) (2'^_2), 



" q-2 ^ q-2 

A2 = ^^Er^(e — l)ln(2 — q). Analytic expressions are correct to leading order in e or to the 
given order in e — 1. The numerical eigenvalues for the fixed points F,G,H are correct up to terms 
0{q-lf. 



21 



Figure captions 



Figure 1: Possible ground states of tlie infinite branched polymer in (a) (b) the contact-rich 
region of the compact phase, (c) the bond-rich region of the compact phase (d) the extended 
phase. 

Figure 2: Numerical phase diagram of the infinite branched polymer in two dimensions in 

the variables (a) jSi — (32 (b) y — t which are defined in the text. The extended phase is 

indicated by E and the compact phase by C. Results are from cluster enumeration (circles) 

PJ33[] and the transfer matrix (open squares) |11| . The full circles give the extended-compact 

collapse transition while the open circles give the compact-compact transition predicted in 

PJ33||. The black squares give the exactly known transition at the percolation point 

and the collapse transition in the strong embedding limit PT|. 

Figure 3: Bethe lattice with coordination number 7 = 3 and central spin ctq. 

Figure 4: p-dependence of the three solutions A (dotted), i?_ (dashed) and C_ (full) of the 

Bethe lattice recursion relations, for y = 2,r = 1.61949. 

Figure 5: Phase diagram of the extended Potts model on the Bethe lattice. The phase 
labelled A corresponds to the extend lattice animal while C corresponds to the compact 
lattice animal. The line indicated by short dashes is first order and the full line is second 
order. The black square at ?/ = 8, r = 2 indicates their meeting point. The jump Ap as 
a function of y between the phases A and C along the first order line is indicated by the 
long-dashed line and is shown on the same scale. 

Figure 6: Fixed point structure which follows from the Migdal-Kadanoff renormalization 
group for the extended Potts model (|2.2|) for g ^ 1 and e ^ 0. The plane p = 1 is also 
shown to guide the eye. 

Figure 7: Phase diagram for the directed branched polymer. The extended phase is indi- 
cated by E and the compact phase by C. 
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